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Abstract. 

The problem of self-trapping of a Bose-Einstein condensate (BEC) and a binary 
BEC in an optical lattice (OL) and double well (DW) is studied using the mean-field 
Gross-Pitaevskii equation. For both DW and OL, permanent self-trapping occurs in a 
window of the repulsive nonlinearity g of the GP equation: g c \ < g < g C 2- In case of OL, 
the critical nonlinearities g c \ and g C 2 correspond to a window of chemical potentials 
/i c i < jj, < fi C 2 defining the band gap(s) of the periodic OL. The permanent self- 
trapped BEC in an OL usually represents a breathing oscillation of a stable stationary 
gap soliton. The permanent self-trapped BEC in a DW, on the other hand, is a 
dynamically stabilized state without any stationary counterpart. For a binary BEC 
with intraspecics nonlinearities outside this window of nonlinearity, a permanent self 
trapping can be induced by tuning the interspecies interaction such that the effective 
nonlinearities of the components fall in the above window. 



PACS numbers: 03.75.Lm, 03.75Kk, 03.75.Mn 
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1. Introduction 

After the experimental observation of a Bose-Einstein condensate (BEC), it is realized 
that a quasi one-dimensional (ID) cigar-shaped trap [Tj is convenient to study many 
novel and complex phenomena. In a cigar-shaped BEC, apart from a harmonic trap, 
double well (DW) [2], periodic optical-lattice (OL) [3], and quasi-periodic bichromatic 
OL [4] traps have been used. Usually, a repulsive BEC is localized in laboratory 
in an infinite trap. However, several other types of counter-intuitive localization of a 
repulsive BEC in a trap of finite height, where possible Josephson tunneling [2], [5] |6] of 
quantum fluids through barriers of finite height is expected to lead to derealization, have 
lately drawn much attention. Of these, self trapping or localization of a repulsive BEC 
predominantly in one site of a DW EJ [10J ED] , triple well [12] and OL [2j [13j HU [15] 
potential has been the subject matter of many investigations. Macroscopic self trapping 
of a BEC was first predicted theoretically [61 [TJ [16] and then observed experimentally 
[21 E]- It is generally believed that self trapping is an intrinsic dynamical phenomenon, 
without any stationary counterpart. In this connection we note that the Anderson 
localization [H [17] in a quasi-periodic or random potential and a gap [H] soliton in a 
periodic potential are both stationary states of the system. In the present critical study 
of self trapping in DW and OL potentials we find that, although, in the former case, is 
is a dynamical phenomenon, in the latter case, contrary to general belief, localization 
takes place in a stationary gap soliton state. There is no symmetry broken stationary 
state corresponding to the self-trapped state in DW. 

To understand the self trapping in an OL and a DW potential, we perform extensive 
numerical simulation of self trapping of a BEC using the solution of the Gross-Pitaevskii 
(GP) equation. In our study, we find a striking similarity between self trappings in 
OL and DW. Both occur in a window of repulsive nonlinearity g of the GP equation: 
g c \ < g < g C 2- The numerical values of g c i and g C 2 depend on the respective trap 
parameters. (The existence of the upper limit of nonlinearity g C 2 was never noted in 
previous studies.) 

Although, the self trapping of a repulsive BEC in a DW represents a relatively 
simple mathematical problem, well described by the analytic two-mode model [71 ITT], 
the self trapping of a repulsive BEC in a periodic OL involving an infinite number of 
wells, on the other hand, poses a formidable mathematical problem and the fundamental 
mechanism in this case is not well understood. For a repulsive BEC in an OL, one has 
localized gap-soliton states [3j [18j [19j [20] . We demonstrate that, in an OL, self trapping 
represents a permanent breathing oscillation of a stable stationary gap soliton. For self 
trapping in OL, we consider compact state(s) localized mostly on a single site of OL, 
in contrast to previous considerations of self trapping on multiple OL sites [13], HU |2T] . 
In the case of an OL, the critical nonlinearities g c \ and g C 2 for self trapping define the 
window of chemical potentials corresponding to the band gap(s). 

We also consider the self trapping of a binary BEC in OL and DW with tunable 
interspecies interaction near a Feshbach resonance [22] . For zero interspecies interaction, 
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self trapping takes place if the nonlinearities g^i = 1,2 of the components are in the 
window g c i < gi < g c2 , and there is no self trapping if g^s are outside this window. In 
the presence of interspecies interaction, the effective nonlinearity of both components 
(arising from a combination of inter and intraspecies interactions) should lie in the 
above window for self trapping. When both g\, g 2 > g c2 , for zero interspecies interaction 
there is no self trapping. We illustrate, using the solution of the GP equation, that, 
for both OL and DW, if g± 2 has an adequate attractive (negative) value, the effective 
nonlinearity of each component is reduced and one could have permanent self trapping. 
When both g\,g 2 < g C 2i f° r zero interspecies interaction there is no self trapping. In this 
case, an appropriate repulsive (positive) interspecies nonlinearity can lead to effective 
nonlinearities of the components in the above window, and hence result in permanent 
self trapping. 



2. Analytical formulation for Double- Well (DW) potential 

We consider a binary cigar-shaped BEC in a DW. Atoms of both species (in two different 
hyperfme states) are assumed to have mass m and number N and the DW acts in the 
axial x direction. Starting with the system of coupled 3D GP equations of the binary 
BEC one can reduce them to the following (dimensionless) ID equations [23] : 

_ [ ^ (X 2 t)] - +^|0,(x,t)| 2 0,(x,t) 
+ g 12 \^(x,t)\ 2 fr(x,t) + V(x)&(x,t), (1) 

where % ^ j — 1,2 denote the species, and wave functions <fii are normalized as 
J^oc \4>i(x, t)\ 2 dx = 1. The suffix x denotes space derivatives and overhead dot time 
derivatives. In ([T]), time t, space x, nonlinearities gi and g± 2 are related to the physical 
observables by [27]: t = u x t phys ,x = x/l, {gi,g\ 2 } = 2{a, ai 2 }\/l with I = ^Ti/mu x and 
A = ujp/uj x is the trap aspect ratio with u x and u p axial and radial (p) frequencies and 
where and a\ 2 are intraspecies and interspecies scattering lengths. The DW is taken 
as [9] 

V(x) =x 2 /2 + Ae~ Kx2 . (2) 
For a single-species BEC, the reduced ID equation is 

i<j>(x,t) = -(l/2)0 x;E (x,t) + (7|0(x,t)| 2 0(a;,t) + V(x)0(x,t). (3) 

In the two-mode model [7J [11] a single- channel BEC wave function <p(x, t) of (J3]) is 
decomposed as [7J 

<j>(x,t) = Mt)®i(x)+Mt)®2(x), (4) 

where spatial modes = 1,2) in the two wells are orthonormalized as 

f$i(x)$j(x) = 5ij and the functions ipi(t) satisfy \ipi(t)\ 2 + \ip2(t)\ 2 = 1. The 
functions ipi(t) are complex and are separated into its real and imaginary parts as 
ipiif) = \i/Ji(t)\ exp(i9i). A population imbalance 

s(t) = (\Mt)\ 2 - l^(t)| 2 )/(|^i(t)l 2 + \Mt)\ 2 ) = (\Mt)\ 2 - \Mt)\ 2 ) (5) 




and phase difference 6 = 6 2 — Q\ then serve as a pair of conjugate variables. The 
approximation (jlj) is then substituted in ([3]), and after some straightforward algebra we 
obtain [7] 



S(t)~- 

6{t) = AS{t) + 2K 



2KyJl-S 2 (t) sin 0(t) 



cos 



where A — gJ dx$>f(x) and 



K = - J dx<&i{x) 



2~ch? 



+ V{x) 



(6) 
(7) 

(8) 



The two-mode equations (jSJ) and (0) are the Hamilton equations [7j S 1 = —dH/d9,9 = 
-dH/dS, for Hamiltonian H = AS 2 (t)/2 - 2K^l - S 2 {t) cos0(t). The transition from 
Josephson oscillation to self trapping happens at H = 2K above a critical A = A c for 

m 



A A c _ 2 + 2 v /l- 1 S 2 (O)cos0(O) 
2K > 2K = 



,S 2 (0) 



(9) 



To perform a numerical calculation using the two-mode model we choose the mode 
functions as [7] 

$i, 2 (x) = [$+(x)±$_(x)]/V2, (10) 
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with the property $i(— x) = $2(^)5 where <&±(x) are the symmetric ground and 
antisymmetric first excited state of ([3]) with potential ([2]) with an appropriate g. The 
parameters of potential (J2J) are taken as A = 16 and k = 10 (used in most of numerical 
calculations below). 

Now we calculate the quantity A/2K for different g and plot in figured] (a). Then 
a critical A c for self trapping is obtained from (Q and plotted in figure [1] (b) versus 5(0) 
for 6(0) = 0. From figure [JJ (b) we see that for 5(0) = 0.1, the critical A C /2K ~ 400. 
From figured] (a) we see that this A/2K is never attained for any nonlinearity g. Hence 
no self trapping can be obtained for 5(0) < 0.1. For 5(0) > 0.2, the critical A C /2K for 
self trapping as obtained from figured] (b), can be attained for g > g\ c as seen in figure 
[JJ(a). However, with further increase of g, A continues to be always greater than A c and 
the self trapping is never destroyed with the increase of g. In our numerical study we 
shall find that the self trapping is destroyed with the increase of g. This is reasonable as 
the two-mode model with the neglect of overlap integrals between the mode functions 
$1,2 (x) is expected to be valid for small g only. A critical nonlinearity for self trapping 
gi c for 6(0) = and different 5(0) is then obtained from the results of A/2K — g and 
A C /2K — 5(0) plots in figures [TJ (a) and (b) and is also plotted in figure [TJ (b). 

3. Numerical results for Double- Well (DW) potential 

The GP equations are solved numerically using the Crank-Nicolson scheme (27] [28] 
with space and time steps 0.025 and 0.0002, respectively, by real-time propagation with 
FORTRAN programs provided in [27] . 

3.1. Single- Channel BEC 

Here we consider DW (]2]) with A = 16 and k = 10, as in some previous studies [9] [10] on 
self trapping. In numerical simulation we create an initial state with a fixed population 
imbalance 5(0), by looking for the ground state at x = Xq of the following asymmetric 
well [9] 

U(x) = (x- x ) 2 /2 + Ae~ KX \ (11) 

The initial 5(0) is achieved by varying the parameter x in ( ITT]) [9]. Once this initial 
state is created, xo is reduced to zero so that DW ( TTTT) reduces to DW (|2]) with xo = 0. 
If Xq is reduced slowly during time evolution, the initial asymmetric state relaxes to 
a symmetric stationary state of DW ([2]) and there is no self trapping. On the other 
hand if Xq is reduced quickly, there is always self trapping provided that g lies in the 
appropriate window: g lc < g < g 2c . This is shown in figure [2] where we plot the density 
at different t after xq is reduced to zero during an interval At = 5, and 20. From figure 
[2] we find that there is self trapping for At = 5 and no self trapping for At = 20, when 
the asymmetric initial state relaxes to a final symmetric state. By actual substitution 
of the wave function of the self-trapped state in the time-independent GP equation we 
verify that there is no stationary state of the DW with this density. The present trapped 
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Figure 2. Density \<j)(x,t)\ 2 versus x at different t when an initiai state, created as 
the ground state of the asymmetric DW (fTTj) . is allowed to evolve in time while the 
asymmetric DW (fTTj) is changed to the symmetric DW ((2]) by reducing xo to zero slowly 
in an interval of time At = 5, 20. When the transition is quick (At = 5) a permanent 
self trapping with an asymmetric profile emerges. For a slow transition (At = 20), the 
system smoothly moves to the symmetric ground state of the symmetric DW ([2]) and 
there is no self trapping. 

wave function was generated from the ground state of DW (II ip and is always positive 
with no node and the present self trapping is dynamical without any static counterpart. 

The self trapping is sensitive to 5(0). There is no self trapping for a small 
5(0)(< 0.1) (viz., two-mode model of section [2]). In figure [3] (a) (upper panel) we plot 
S(t) versus t to illustrate the trapping and oscillation for different g (= 1, 10, 100, 1000) 
for 5(0) = 0.3. For g — 1, there is Josephson oscillation with (5(t)) = 0. For g = 10 
we have (S(t)) ~ 0.3 insuring self trapping. As g is further increased, the nonlinear 
term in the GP equation becomes much larger than the Gaussian wall Aexp(— kx 2 ) 
in ([2]). Consequently, the DW appears like a single well and classical center-of-mass 
oscillation appears for a large g. For g = 100 the oscillation is irregular because the 
Gaussian wall is not fully negligible. Finally, for g = 1000 the Gaussian wall can be 
completely neglected and we have clean classical center-of-mass oscillation. The S(t) 
versus t plot for 5(0) = 0.3 and g — 10 at large times is shown in the lower panel of 
figure [3] (a), where self trapping is confirmed up to t = 5000 — an interval much larger 
than tunneling time which is at best ~ 100. The evolution of S(t) at large times 5(oo) 
for different 5(0) as g is increased is shown in figure |3](b). A nonzero 5(oo) ensures self 
trapping, which appears in the window: g c \ < g < g c2 . The prediction of the analytical 
two-mode model of section [2j from figure [1] (b), for g c \ are shown by arrows in figure 
|3] (b) for 5(0) = 0.3 and 0.6 in good agreement with numerical simulation. It is 
noted that in all calculations reported here g C 2 < 100 — a medium value of nonlinearity 
where the mean-field GP description is well justified. For larger values of nonlinearities, 
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Figure 3. (a) (upper panel) Population imbalance S(t) versus t for 5(0) = 0.3 in 
DW ([2]). (lower panel) Long-time S(t) versus t for g = 10. (b) Average population 
imbalance at large time 5(oo) versus g for 5(0) = 0.3 and 0.6. The arrows are results 
for g c i from figure [T] (b) . 

beyond mean- field corrections to the GP equation [29] could be relevant [30]. However, 
we expect the general conclusions of this study to remain valid for the present set of 
parameters. 
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Figure 4. (a) (upper panel) Population imbalance S(t) versus t for a binary BEC with 
5(0) = 0.43 and nonlincarities g x = g 2 = 50 and g 12 = 0, -40 and -49.9 with DW ©. 
(lower panel) Long-time S(t) versus t dynamics for 512 = —40. (b) Average population 
imbalance S(oo) at large time versus effective nonlinearity g\ + 512 for S*(0) = 0.43. 
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3.2. Binary BEC 

We consider the self trapping of a binary BEC in a DW for g$ > g C 2 and gi < g c \. 
In both cases no self trapping is possible for zero interspecies interaction (g\ 2 = 0). 
Of these, for g^ < g c \ we need a repulsive (positive) g 12 to make both components 
sufficiently repulsive to have self trapping. But this case usually leads to a uninteresting 
stationary phase separated configuration [31] with the first species occupying the first 
well of the DW and the second species occupying the second well. A more interesting 
situation emerges for g$ > g c2 , when we require an attractive (negative) gi2 to make 
both components appropriately repulsive to have self trapping. However, self trapping 
appears only for an intermediate S(0) for an appropriate gi and g 12 . To illustrate, we 
consider the symmetric case with gi = g 2 = 50 in ([T]) and (T5]) and consider an attractive 
(negative) gi 2 . In the binary case, the self trapping is not so good for k = 10, while 
the barrier between the two wells in (T5]) is narrow and we take k = 5 in the numerical 
simulation. Permanent self trapping is found to occur for a narrow window of 5(0) 
around 5(0) = 0.4. For S(0) > 0.5 and S(0) < 0.3 self trapping for a small interval of 
time could be obtained. 

The initial state with S(0) = 0.43 for both components is created using DW ( TTTT) 
and reducing x to zero. In the upper panel of figure [4] (a), we plot S(t) versus t for 
both components for three values of g± 2 . For g± 2 > 0, no self trapping is possible and 
we are in the domain of irregular oscillation as shown for g 12 = 0. For an attractive gi 2 , 
there is permanent self trapping as illustrated for gi 2 = —40. With further increase of 
|<7i2| the self trapping disappears and regular sinusoidal Josephson oscillation appears 
as shown for g 12 = —49.9 in figure |4] (a). The lower panel of figure H] (a) illustrates 
the population imbalance S(t) for g\ 2 = —40 at large times, which confirms robust self 
trapping. In figure @] (b) we plot S(oo) for both components versus g\ + g\ 2 . In the 
symmetric case with g\ = g 2 , the wave function of the two components are equal and 
the quantity g x + (712 is the effective nonlinearity of each component. Figure IU(b) shows 
that self trapping occurs for the window of effective nonlinearity 0.2 < g\ + g\ 2 < 20 
for S(0) = 0.43. For any given initial population imbalance and for either sufficiently 
small or sufficiently large effective nonlinear interaction strength gi + g 12 , the system 
is in the oscillation regime. For small values of effective interaction one has Josephson 
oscillation and for large values of effective interaction one has the classical center-of- 
mass oscillation. For intermediate interaction strength, the system may make transition 
to self trapping for an appropriate S(0). 

4. Analytical formulation for Optical-Lattice (OL) potential 

For an OL, Eqs. ([!]) and (j3]) remain valid but with the variables bearing the following 
relations to the corresponding physical observables [24]: t = (7r/L) 2 (^/m)t p hy S , x = 
xn/L,Vo = m(L/iih) 2 V i, {gi, gi 2 } = (2NLmu±/nh){ai,ai 2 }, where L the wavelength 
and Vq the strength of OL potential: Vql(x) = — Von 2 h 2 / (mL 2 ) cos(2nx/L). The 
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Figure 5. (a) Chemical potential fj, versus Vo showing the bands (shaded area) and 
gaps (white area) for OL (TT^j) . The hatched area between the lower line with crosses 
and the upper line with pluses represents the variational domain of the appearance of 
Gaussian gap solitons. The full white area above the lower line with crosses denote 
the exact domain of the appearance of gap solitons. (b) The critical values g c \ and g C 2 
defining the window of nonlinearities of Gaussian gap solitons versus Vq. 



dimensionless OL potential is given by 

V(x) = -F cos(2x). (12) 

To critically examine the self trapping in a OL where the BEC is localized in one of 
the sites of the OL, we present a Gaussian variational analysis. The stationary state is 
described by ([3]) with the time derivative term i<p replaced by [Mp where /i is the chemical 
potential. The Lagrangian for that stationary equation is given by [19] 

/i0 2 (x) - (f) 2 x (x)/2 - g(j) 4 (x)/2 + Vocos(2x)(j) 2 (x) dx - //. (13) 

-oo L - 1 

This Lagrangian can be analytically evaluated by considering a simple form for the wave 



function 0. Using the Gaussian form [25] (j)(x) = it 1 / 4 y7V/u> exp[— x 2 /(2w 2 )] with J\f 
the norm and w the width, the Lagrangian can be written as 

L = ~ 1) - f~ 2 + Vo^eM-™ 2 ) ~ t^-- (14) 
The variational equations [2S] dL/dfx = dL/dw = dL/dM = yield M — 1 and 

1 1 9W -4Vow 4 exp(-u> 2 ), (15) 



V2tt 

^ = A + ~7f- ~ K,exp(- U ; 2 ), (16) 



which determine the width and the chemical potential. We have set M = 1 in ( TT5"j) and 
( TT6|) after derivation. 

It is instructive to study the band and gap structure of the Schrodinger equation 
in the periodic OL potential (Jl2|) : (J,<fi(x) = —(l/2)<p xx (x) — V cos(2x)(j)(x) [26]. In a 
periodic potential the quantum excitation spectrum of a system consists of bands and 
gaps. The bands allow unlocalized plane wave solution modulated by a periodic function 
with the same period as the periodic potential known as Bloch wave. The gaps permit 
localized solution [26]. However, the single-particle linear Schrodinger equation does 



Self-trapping of a binary Bose-Einstein condensate induced by interspecies interactionll 
0.8 i L 




X 



Figure 6. Density \cj)(x,t)\ 2 versus x at different t when an initiai Gaussian state 
with trapping measure D(0) = 0.99 is released at x = t = on the OL. The variational 
result for the corresponding gap soliton is shown by crosses. The positive part of the 
potential is also shown in arbitrary units. 

not allow any solution in the gap. The band (shaded regions) and gap (white regions 
above the lowest shaded region) of the spectrum of the Schrodinger equation with OL 
potential ([12]) is shown in figure [5] (a). The nonlinear ID GP equation fl3]) for repulsive 
interaction (positive g) permits localized solutions in the gaps, called gap solitons, where 
the chemical potential /i lies in the gap. In the gaps, localized gap solitons are possible 
in the presence of an appropriate nonlinearity g. 

We now find the condition for a gap soliton for a positive (repulsive) g by directly 
solving ( !T5|) and (fl6|) for all g and V . As fi tends to the upper edge of the lowest band, 
one has the lower limit g c \ of the formation of a gap soliton, denoted by a line with 
crosses in figure [5] (a). As the repulsive nonlinearity g is increased, ( I15jl and ( Tl6l) permit 
solution up to a maximum value g& which determine the upper limit of the formation 
of a gap soliton denoted by the line with pluses in figure [5] (a). The area between the 
two lines determine the domain in which the Gaussian gap solutions are allowed. (Gap 
solitons of non-Gaussian shape, possibly occupying many OL sites, are possible in the 
whole white region above the lowest band gap in figure |5] (a).) The nonlinearities 
corresponding to these two lines can be calculated using ( 115]) and ffT6]) and yields the 
critical g c \ and g C 2 for self trapping. These nonlinearities are plotted in figure [5] (b). 

5. Numerical results for Optical-Lattice (OL) potential 

5.1. Single- Channel BEC 

The numerical simulation is started by releasing a Gaussian state in the center of the 
OL during time evolution of the GP equation. The width of the Gaussian is taken such 
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Figure 7. (a) Trapping measure D(t) versus time t for Vq = 5 for different 5. 
(b) Trapping measure at large time D(oo) versus g for different Vq — 1,2,3,5. The 
respective arrows represent nonlinearities g c \ and g c i from figure [5] (b). 



that the initial state stays mostly in the central well of the OL. An estimate of self 
trapping is given by the following function called trapping measure 

D(t) = I"' 2 \<p{x,t)\ 2 dx. (17) 

J-ir/2 

Here we note that it is the wave length of the OL — Vq cos(2x). Hence the function D(t) 
determines the matter inside a single OL site. In case of ideal self trapping in a single 
site D(t) = 1, and D(t) will tend to zero when self trapping is fully destroyed. 

Self trapping occurs easily if the nonlinearity g is appropriate for a Gaussian gap 
soliton (viz. figure |5]). We take an initial Gaussian state with D(0) = 0.99 and release 
it on the OL with Vq = 3 and g = 5. The density profile of the trapped state are shown 
in figure EJ The small dispersion of the density profile at different t guarantees good 
self trapping. We also plot the density of the variational gap soliton, in good agreement 
with the self-trapped state, in figure [6] 

In Fig [7] (a) we plot D(t) versus t for Vq — 5 and different g. For g < 10, D(t) 
remains close to 1 for t < 100. However, if we continue to large t(< 10000) (a time 
much larger than the tunneling time of few hundred as in Fig [7] (a)), D(t) remains close 
to 1 for a window of nonlinearity 0.01 < g < 10 denoting permanent self trapping. For 
illustration, in figure [7J (b), we plot the trapping measure D(t) at large times D(oo) 
versus g, where D(oo) is non-zero in the window g c \ < 9 < 9c2 corresponding to 
permanent self trapping and is zero outside showing no self trapping. In figure [7J (b) 
there are two arrows for each Vq corresponding to g c \ and g C 2 as obtained in figure [5] (b). 
The domain between the two arrows representing the region of allowed gap solitons (see, 
figure [5]) agrees well with the domain of self trapping represented by non-zero D(oo) 
as seen in figure [7J (b). From this fact and also from the proximity of the variational 
solution for density of a gap soliton with that for the numerical self-trapped state in 
figure El we conclude that the self trapped state represents small oscillation of a stable 
gap soliton. In this connection, the finding in Ref. [13], that a self-trapped state in an 
OL is always temporary, is not fully to the point. They missed the fact that for the 
window of nonlinearity g c i < 9 < 9c2 the self trapping could occur in a stable stationary 
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Figure 8. (a) Trapping measure D(i) of a binary BEC versus time t for different 
<?i — 5i2- The curves are labeled by g\ — gn and the component number in parenthesis, 
(b) Trapping measure at large time -D(oo) versus g\ — gn for different Vq. The arrows 
are g c \ and g c i from figure O (b) . 



gap soliton state leading to a permanent trapping. Just outside this small window of 
nonlinearity a temporary self trapping at small times may take place, as can be seen in 
figure [7| (a). 

5.2. Binary BEC 

To study self trapping of a binary BEC we consider large repulsive intraspecies 
nonlinearities g\ — gi — 50, which do not allow self trapping for zero interspecies 
interaction g\ 2 = for Vq = 1,2,3,5, as seen in figures El (a) and (b). If we introduce 
an attractive (negative) interspecies interaction g 12 , then in each channel the effective 
nonlinearity will be reduced and for a sufficiently large and attractive gi2 one can have 
self trapping. The trapping dynamics for this system is illustrated in figure M (a) where 
we plot trapping measure D(t) versus t for different g 12 and V — 5,g± — g^ — 50. The 
results for the two components are practically the same in most cases. The self trapping 
appears for a small (gi + gi 2 ) (~ 0.01) and disappears for large (gi + gi 2 )(> 10). In this 
symmetric binary BEC, (g± + g\ 2 ) provides a good measure of the effective nonlinearity 
controlling self trapping. In figure M (b) we plot the trapping measure at large times 
.D(oo) versus effective nonlinearity g + g± 2 of the binary BEC for Vq = 1,2,3,5. The 
plots of figures El (b) and El (b) are qualitatively quite similar, showing that (gi + gi 2 ) 
is a good measure of effective nonlinearity of the binary BEC. However, if |g 12 | is taken 
to be larger than g then the effective nonlinearity becomes attractive corresponding to 
a negative g\ + g\ 2 . This domain of nonlinearity corresponds to permanent symbiotic 
bright soliton [32] and consideration of self trapping is inappropriate. 

Finally, we consider a binary BEC with small intraspecies nonlinearities and zero 
interspecies nonlinearity (g 12 = 0), that does not allow self trapping, as found in figures 
El (a) and (b). As an illustration we consider the zero values for the intraspecies 
nonlinearities: g\ = g 2 = 0. In the presence of appropriate repulsive interspecies 
nonlinearity (g\ 2 > 0), one can have self trapping as illustrated in figures [9] (a) and 
(b). In figure [9] (a) we plot the trapping measure D(t) versus t for different repulsive 
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Figure 9. (a) Trapping measure D(t) of a binary BEC versus time t for Vb = 5 and 
different g\2- Curves are labeled by g\2 and component number, (b) Trapping measure 
at large time D(po) versus 512 for different Vq. The arrows represent nonlinearities g c \ 
and g C 2 from figure [5] (b). 

(positive) interspecies nonlinearity g 12 and V = 5, gi = g 2 = 0. In this case g\ 2 plays 
the role of effective interaction. There is a window of g\ 2 values g c \ < g\ 2 < g c2 with 
g c i ~ 0.01 and g c2 ~ 10, where permanent self trapping can be achieved. In figure [9] (b) 
we plot large-time trapping measure D(oo) versus g± 2 for different Vq and g\ — gi = 0. 
Qualitatively, this plot is quite similar to those in figures [7] (b) and M (b) showing the 
universal nature of these plots. 

6. Summary and Discussion 

We demonstrated that self trapping of a BEC or a binary BEC without interspecies 
interaction in OL and DW occurs for a window of repulsive intraspecies nonlinearity g 
(Oci < 9 < 9c2), where g c \ and g c2 depends on the trap parameters. For a binary BEC 
with the intraspecies nonlinearities outside this window, a self trapping can be induced 
by a non-zero interspecies nonlinearity gi 2 such that the effective nonlinearities fall in 
this window. For intraspecies nonlinearities gi below g c \ (g^ < g c i), this is achieved 
by an appropriate repulsive (positive) gi 2 . For intraspecies nonlinearities gi above g c2 
(9i > 9c2)i one can have self trapping by introducing an appropriate attractive (negative) 
gi 2 . In case of self trapping in an OL, the permanently self trapped state represents 
breathing oscillation of a stable stationary gap soliton. On the other hand, the self 
trapping in a DW is purely a dynamical phenomenon without any underlying stationary 
state. However, the self trapping of a BEC and a binary BEC in both DW and OL could 
be permanent. 

In previous studies of self trapping of a BEC, the existence of a lower limit g c \ of 
nonlinearity was noted [2j [7J [13j [H]. However, the disappearance of self trapping above 
an upper limit was not realized. (A similar upper limit appeared in a study of a Fermi 
superfluid in a DW |9j.) In a previous study of self trapping in an OL, in contradiction 
with the present investigation, it was concluded that self trapping was always transient 
and should disappear at large t [T3j . 

The principal findings of this critical study of self trapping are the following, (a) 
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The self-trapped states in an OL potential are essentially the stationary gap solitons. 
The self trapping in DW potential is entirely dynamical in nature and there are no 
stationary states in this case, such as the gap solitons of the OL potential. (There is 
no periodic potential and no band and band gap, specially for the case of a DW with a 
shallow barrier, which we shall study here.) It was generally believed that self trapping 
in both OL and DW potentials are dynamical in nature, (b) Self trapping is stopped 
beyond an upper limit of interaction in both cases, (c) In case of coupled systems 
with inter-species interaction self trapping is possible for domain of intra-species atomic 
interaction where no self trapping is allowed in uncoupled systems. 

In the experiment of self trapping in an OL and in related theoretical studies a 
localized state over tens of OL sites [131 EH] was considered in contrast to that in 
predominantly a single OL site. Such states extended over multiple OL sites could 
possibly be a combination of multiple compact gap solitons. In another study, such 
states have been suggested to be a new type of spatially extended state in the gap [2T] . 
The compact self trapped state on OL considered in this paper are different from the 
spatially extended states considered in other studies. Nevertheless, with the present 
experimental control over a BEC, it would be possible to study self trapping of the 
compact states on OL as considered in this paper. 
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